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ABSTRACT 

Context. Most observational studies so far point towards brown dwarfs sharing a similar formation mechanism as the one accepted 
for low mass stars. However, larger databases and more systematic studies are needed before strong conclusions can be reached. 
Aims. In this second paper of a series devoted to the study of the spectroscopic properties of the members of the Lambda Orionis Star 
Forming Region, we study accretion, activity and rotation for a wide set of spectroscopically confirmed members of the central star 
cluster Collinder 69 to draw analogies and/or differences between the brown dwarf and stellar populations of this cluster. Moreover, we 
present comparisons with other star forming regions of similar and different ages to address environmental effects on our conclusions. 
Methods. We study prominent photospheric lines to derive rotational velocities and emission lines to distinguish between accretion 
processes and chromospheric activity. In addition, we include information about disk presence and X-ray emission. 
Results. We report very large differences in the disk fractions of low mass stars and brown dwarfs (~58%) when compared to higher 
mass stars (26^ %) with 0.6 M Q being the critical mass we find for this dichotomy. As a byproduct, we address the implications of 
the spatial distribution of disk and diskless members in the formation scenario of the cluster itself. We have used the Ha emission 
to discriminate among accreting and non-accreting sources finding that 38^ % of sources harboring disks undergo active accretion 
and that his percentage stays similar in the substellar regime. For those sources we have estimated accretion rates. Finally, regarding 
rotational velocities, we find a high dispersion in vsin(i') which is even larger among the diskless population. 

Key words. Stars: formation - Star: low-mass, brown dwarfs - open clusters and associations: individual; Collinder 69 



1. Introduction 

This is the second paper of a series devoted to studying, from a 
spectroscopic point of view, the young population present in sev- 
eral associations belonging to the Lambda Orionis Star Forming 
Region (LOSFR). In this paper, we systematically analyze sev- 
eral properties of the confirmed members of the central cluster 
Collinder 69 (C69); including the presence of disks, accretion 
onto the central star/brown dwarf, rotation and activity. 

(from now on Paper I), we presented 



In Bayo et al. (2011 



a very detailed and complete spectroscopic census of Collinder 
69; the oldest (~ 5 — 12 Myr) of the associations belonging to 
the LOSFR. In short, this star forming region is located at ~400 
pc ( Murdin & Penston|1977| l representing the head of the Orion 
giant. Its center is dominated by the 08III multiple star A Orionis 
and comprises both recently formed stars from 0.2 M to 24 M Q 
and dark clouds actively forming stars. 

The main goals of this paper are to study in detail the prop- 
erties of the C69 confirmed members (e.g. rotation, activity, disk 
fractions and accretion rates), analyze the different populations 



* Based on the ESO observing programs 080.C-0592 and 078.C- 
0124; and observing programs from Calar Alto, Keck, Subaru and 
Magellan. 



within the cluster, and compare our results with other low mass 
star forming regions of similar and different ages. In fact, the 
evolutionary status of C69 seems to be specially suited to the 
study of disk evolution and accretion at the low end of the mass 
spectrum. According to our observations, in almost every bin 
in mass (for masses lower or equal to 0.6M o , ~3700K assum- 
ing a 5 Myr isochrone from |Baraffe et al. 1998| l, we find a 



diversity of disk harboring sources: from those with optically 
thick disks undergoing active accretion (onto the central star) to 
others that have completely lost their primordial circumstellar 
disks. Finally, and as a byproduct of this study, we try to put 
our findings in context of the current theory of how the LOSFR 
as a whole was formed (triggered by a supernova, see Dolan & 
Math ieu|2001| and more discussion on this scenario on Paper I). 



This paper is organized as follows: In Section [2] we provide 
a description of the data analyzed. In Sections [3] to |5] proper- 
ties such as rotation velocity, activity levels, accretion processes, 
disks, variability and spatial distribution of the population of 
C69 are studied. And, finally, our summary conclusions are pre- 
sented in Section [6] 

Unlike in Paper I, we have grouped the most interest- 
ing/puzzling objects into different categories (following the sub- 
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sections of Sections [3] and [4]) and discuss their peculiarities in 
Appendix A. 



2. Observations and data analysis 

In this work we make use of an extensive spectral database 
that our group has been gathering during more than seven years 
and that is described in Paper I. Furthermore, we make use of 
the measurements and parameters derived for confirmed mem- 
bers published in |Dolan & Mathieu| (|T999] |200l); |Barrado y 
|Navascues et al.| (|2007|>; ISacco et al.| (|2008|); |Maxted et al. 
d 2008b ; |Morales-Calder6n|(|2008|l; |Bayo| ( |2009| |; |Barrado et aT 
pOTI) and |Franciosini & Sacco| ( |2011 1 ). 



2.3. Fundamental parameters database 



2.1. Spectroscopic database 

In total, we analyzed spectra obtained by us corresponding to 
~100 confirmed members, with several objects observed more 
than once, plus data corresponding to ~70 objects that were stud- 
ied in at least one of the papers cited before and for which we do 
not have our own observations. 

In short, our own database comprises optical spectra of 
confirmed members of C69 with temperatures between ~2800 
and ~4700 K. The spectra have resolutions in the range 600- 
1 1250, wavelength coverages from ~5000 A up to ~10400 A and 
have been obtained with different instruments and telescopes at 
Mauna Kea, Las Campanas, CAHA and VLT. 

The data reduction of most of our campaigns was performed 
with IRAfQ following the standard steps. For the analysis of the 
spectra, motivated by the large amount of data, we developed a 
tool that provides for any given set of lines in an automatic man- 
ner (instrumental response corrected) full width half maximum 
(FWHM), full width at ten per cent of the flux (FWi %), and 
equivalent widths (EWs). The description of this routine and a 
study on the effect of the resolution on parameter determination 
is given in Appendix A of Paper I. 



2.2. Photometric and X-ray database 

In addition to the spectroscopic database, we have made use of 
the photometric dataset described in Paper I and analyzed in 
detail in |Morales-Calder6n| ( |2008] l; |Bayo| p009) and Morales- 
Calderon et al 2012 in prep., and of the X-ray data presented in 
Barrado et al.|(2011 1 and Franciosini & Sacco] ( |201 1\ . The pho- 



tometric database includes photometry from the optical to the 
mid-infrared (MIR) 

Taking a 5 Myr isochrone as a reference ( Baraffe et al.|1998 



Chabrier et al. 2000 Allard et al. 



2003| l, the completeness of 



the optical data is located at ~20 Mj up for the whole cluster and 
even down to ~3 M Jup for a field to the East of the center of the 
cluster. In the near infrared (NIR), the data is complete down to 
~10 Mj U p in almost the whole cluster and down to ~30 Mj up in 
the outskirts. For the MIR, our database is complete at 3.6 /im 
down to -40 M Jup . 

Finally the X-ray data is complete down to ~O.3M0 with 
detections in Barr ado et al.| ( |201 l| l for confirmed members as 
low as ~0.1M G . 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc. under contract to the National Science 
Foundation. 



We have adopted the results from Morales-Calderon ( 2008 1 and 



Barrado y Navascues et al. ( |2007 1 regarding the presence of cir- 
cumstellar disks and substellar analogs. Objects with NIR and/or 
MIR slope above photospheric values are considered to harbor 



disks. This slope is analyzed as in Lada et al. (2006) based on 
the a parameter. 

The photometric database described before was used to feed 
VOSA ( |Bayo et al.||2008| ) and determine, via Spectral Energy 
Distribution (SED) fit, the effective temperatures and the bolo- 
metric luminosities for all confirmed members. These two pa- 
rameters were used in Paper I, along with theoretical isochrones 
from the Lyon gro up ( |Baraffe et al.|1998[ |Chabrier et al.|2000[ 
|Allard et al.|2003| l, to determine masses for the confirmed mem- 
bers. 

For the mass determination we assumed an age of 5Myr for 
C69 since that was the best fit to the produced HR diagram. 
We note, however, that in Paper I we also set an upper-limit for 
the age of C69 of 20Myr. We will address this point again in 
Section [ 



3. Rotational velocities 

For the sample of sources observed with the highest resolution 
(those observed with Magellan/MIKE, R~11250; 14 sources in 
total), we estimated their projected rotational velocity. We add to 
this sample the objects with v sin(/) determination in |Sacco et al. 
(2008) to study the rotation of ~25 members of C69. 

3.1. Technique 

We based our determination of the projected rotational velocity 
on the comparison of the observed spectra with Kurucz models 



( Caste 111 et al.||1997| ) synthesized for different rotational veloc- 
ities. Since our estimation of v sin(j) is based on comparisons 
with models of a specific T e ff, and there is a known dispersion 
between effective temperatures based on spectral types and those 
derived using models (T e g scales, see Paper I), we followed a 
three step process: 

1. First, we derived T e ff and log(g) by comparison with the- 
oretical models. We built a grid of synthetic spectra using 
different collections: Kurucz (Castelli et al. 1997), NextGen 
( jHauschildt et al.|1999| l and Dusty ( | Allard et al. 12003] ) cover- 
ing effective temperatures in the range 5000-2000 K, assum- 
ing solar metallicities, and a range of log(g) between 3.5 and 
5.0. We lowered the resolution of both our spectra and the 
Kurucz models to R~200, which is the resolution provided 
in the public servers by the Lyon group, in order to perform 
a direct fitting process. Table [T] shows the best fitting set of 
parameters for each case. 

2. Next, we prepared a second grid of high resolution Kurucz 
models (the ones we could generate with the same resolu- 
tion as our Magellan/MIKE spectra), for different values of 
vsin(/) for direct comparison. 

For the comparison we chose the wavelength range from 
6090 A to 6130 A because the signal-to-noise ratio (SNR) of 
our spectra was not very high in all cases, and, for the range 
of temperatures to be analyzed, two of the most prominent 
photospheric absorption lines (from Ca I) are located within 
this region. Besides, for the grid of high resolution models, 
we assume log(g) of 4.0 dex that is a suitable value (accord- 
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ing to models) for these young cool confirmed members of 
C69. 

3. Since our high resolution grid of models did not include 
models from the Lyon group, we only analyzed sources 
where Kurucz was the best choice in Step 1, or those for 
which the difference in reduced x 1 between Kurucz and 
NextGen was not significant (LOri050, LOri055, LOri057, 
LOri060 and LOri075). 

We performed a model fit to the Ca I lines described in Step 
2 (see Fig.[T| in order to estimate the rotational velocities of 
our sources. Prior to this, we included the instrumental re- 
sponse measured on the MIKE lamp arc. The mean FWHM 
value measured on the lamp spectra was 0.4914 A, corre- 
sponding to a velocity of ~20 km/s (the lower limit we can 
provide for several cases). 

3.2. Main results 

The estimated rotational velocities from this work are provided 
in Table[T[ a nd also in Table[6]where we add those determined by 
Sacco et aL|(|2008|l. In general, the range of values determined in 



both studies agree very well. In particular, there are three sources 
for which both Sacc o et aX] ([20Q8 ) and us have estimated vsin(/). 
In two of these cases, LOri055 and LOri057, both estimations 
are in very good agreement, but in the case of LOri060 we pro- 
vide a value ~10 km/s higher than the one obtained by Sacco 
et al.| ( |2008] l. For this source, our determination of the rotational 
velocity suffers from the fact that the S/N of the spectrum was 
the lowest of the sample and, for this reason, the wings of the Ca 
I doublet are not so "clean" as in the other cases. Therefore, we 
take their upper-limit of ~ 20 km/s as a more robust measure- 
ment. 

In Fig. |2j we display the projected rotational velocity as 
a function of spectral type for those objects where a value of 
vsin(/) (either calculated by us or from Sacco et al. 2008 1 is 
available. Although a significant fraction of the measurements 
are upper limits, we can see how, as expected given the youth of 
the sample, the general trend for objects in C69 is to rotate faster 
than the old disk population of low-mass stars and brown dwarfs 
from Mohanty & Basri (2003 ) (starred symbols). 

Furthermore, although the dispersion of vsin(/) values 
among the members of C69 is large; if we consider the NIR and 
MIR slope as a proxy for disk presence, we see how this disper- 
sion is larger among the diskless population than among disk- 
harboring sources. Although we are dealing with small num- 
ber statistics, this could be the result of the disk locking effect 
( Bouvier et al.p997) , that explains the dispersion in rotational 
velocities as a result of different decoupling timescales between 
the star and the disk. 



Table 1. T e ff, log(g) and rotational velocities derived by com- 
parison with models for the sample of sources observed with 
Magellan/MIKE. We also display other relevant information 
about the sources such as accretion rate (derived in this work) 
or binarity. 



Source 


Mdl 1 


SpT 


T e ff 


log(g) 


vsin(;') 


A 1 

Acc 


log(M ncc ) Binary J 






(K) 






(km/s) 




(M Q /yr) 


LOri017* 


Kur 




4250 


4.0 


70 






LOH031 


Kur 


M3.5 


3750 


4.0 


40 


N 




LOri042 


Kur 


M4.0 


3750 


4.0 


30 


N 




LOri050 


NG 


M4.0 


3700 


4.0 


60 


Y 


-10.878±0.05 M08, S08 


LOri055 


NG 


M4.0 


3700 


4.0 


<20 


N 




LOri056 


NG 


M4.5 


3700 


4.0 




N 




LOri057 


NG 


M5.5 


3700 


4.0 


<20 


N 




LOri058 


NG 


M3.5 


3500 


4.0 




N 




LOH059 


NG 


M4.5 


3500 


4.0 




N 




LOri060 


NG 


M4.0 


3500 


4.0 


20-40 


N 




LOri063 


NG 


M4.0 


3700 


4.0 




Y 




LOri068 


NG 


M4.5 


3700 


4.0 




N 




LOri073 


NG 


M5.0 


3700 


4.0 




N 




LOri075 


NG 


M5.0 


3400 


4.0 


65 


N 


M08, Bll 



'There were problems on the extraction of several orders. We have 
used the T eff corresponding to the SED fit provided in Paper I. 
1 Model collection corresponding to the best fitting model: NG = 
Nextgen; Kur = Kurucz 



2 Y indicates that the object fulfills the accretion criterion by |Barrado y 
|Navas cues & Martin ( 2003 ) (those marked with N do not) 



s Binary system according to Maxted et al. (2008 , M08), |Sacco et al. 
p008l S08) or this work (Bll). 



L0ri031 




6090 



6100 



6110 
Lambda (A) 



6130 



6130 



Fig. 1. Example of vsin(/) estimation for LOri031. In black we 
have plotted the science spectrum and on top of it with differ- 
ent colors the Kurucz model synthesized for 3750 K, log(g) of 
4.0 and five different velocities. The synthetic spectra have been 
produced with the same resolution as the observed one. 



4. Activity and accretion 

Activity, accretion, and mass loss processes can be studied 
through the analysis of emission lines in the spectrum of young 
stellar and substellar objects. We can group these emission lines 
into two categories: (i) The forbidden emission lines of [OI] 
(,15577, 6300, 6364 A), [SII] (/16717, 6731 A) and [Nil] (,16548, 
6581 A) that have been attributed in the literature to low den- 
sity regions such as winds, and are therefore a tracer of the 
mass loss process (see Shu et al. 1994 Hartmann et al. 1994 
and |Hartmann|| 1 999) . (ii) The permitted emission lines of He 



I 0t6678 A), Ha (i6563 A) and the Call triplet (^8498, 8542, 
8662 A) which are characteristic of classical T Tauri stars and 
their substellar analogs and trace accretion processes, although 
they are also known to be signposts of chromospheric activity 



(Mart in et al.|2001||Natta & Testi|2001||Natta et al.pMl |Testi| 
et al.||2002| |Mohanty et al.||2003|lj ayawardhana 1Fal||2002b|a[ 



2003a b Barrado y Navascu es et al.|2002||2003HWhi te & Basr 
|2003D . 

Given the evolutionary status of the members of C69 and 
the environment where they are located, most of the analyzed 
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Fig. 2. Spectral Type vs vsin(z) for confirmed members of C69 
from this work and Sac co et aL] ( |2008) (black dots). For com- 
parison, we include the old disk population brown dwarfs from 
Mohanty & Basri (2003) with five pointed star symbols. From 
the C69 members we have highlighted those exhibiting peculiar- 
ities: large open squares surround spectroscopic binaries, large 
open circles represent objects with infrared excess (assumed to 
be a signpost of the presence of a disk) and red dots have been 
plotted on top of those objects classified as accretors accord- 
ing to the saturation criterion by Barrado y Navascues & Martin 
(2003 ) (see text for details). Finally LOri075 is highlighted with 
a label and is discussed in Appendix A. I' 



spectra show a rich variety of emission lines. We summarize the 
results from our analysis in the following subsections. 



sources are in the BD domain and the rest have masses above the 
hydrogen burning limit. 

Further comparison of the He I and Ha EWs of the 24 
sources with He I detection shows a clear correlation among 
them; the stronger the Ha emission, the stronger the He I one; 
but our sample is not statistically significant to derive any quan- 
titative relationship between them. 

The fact that we detect He I mainly in stars and not in brown 
dwarfs is not surprising since the typical structure of the He I line 
is much narrower than that of Ha, for example, and therefore this 
emission can be hidden if the spectrum does not have high SNR 
and resolution. 

In addition to this technical caveat to detect He I, the ex- 
tended nebular emission of the cluster has to be taken into ac- 
count when dealing with fiber spectrographs. In particular, in 
two cases from the VLT/FLAMES spectra, the helium line could 
be contaminated by nebular emission (see in Fig.[3]the distribu- 
tion of the sky fibers and in Section 4.3 the quantification of this 
contamination for the Ha line). 



Q 9.8 



9.4 - 
84.4 



A * 

3K O 



SO ft'. 



.a? •• 



o B&C 

* CAFOS 

* FLAMES 
o LRIS 

a MIKE 
□ TWIN 



84.2 84.0 83.8 83.6 83.4 83.2 

R.A 



4.1. Spectroscopic emission lines among the C69 sample 

For our sample, the forbidden emission lines (given their nar- 
row nature) were only detected in some of the higher resolution 
spectra (see tables [2] and [3J, but since most of these detections 
are quite marginal we cannot be certain whether the signal comes 
from the source itself or from the nebular environment. 

The Ca II triplet was covered in very few of our runs and 
only detected in one object (C69-IRAC-005) that is analyzed in 



detail in Appendix A. 2. 2 



On the other hand, we have measurement/s of Ha for 156 
sources out of the 172 spectroscopically confirmed members. 
The vast majority of these sources(140 objects, ~90% of the 
sample) show the line in emission. The EW of this line is used 
to distinguish between chromospheric or accretion origin of the 
emission in Section |4"31 

Regarding He I; we have detected this line in 24 objects, all 
in emission, and we provide the EWs in tables 2]and[3] All these 
sources present also intense Ha emission (with a minimum EW 
of ~5A) and six of them are classified as accretors in Section 4.3 



Regarding the mass range where He I is detected; only six of the 



Fig. 3. Contours corresponding to an Ha image of C69 (ex- 
tracted from the Ha six arcmin resolution all sky survey com- 
piled by Finkbeiner 2003). We have plotted the location of the 
C69 members using different symbols (all of them in black, see 
legend of the figure) according to the instrument with which the 
individual spectrum was obtained. As can be seen a significant 
fraction of the sources have been observed more than once by us. 
The locations of sky fibers in the VLT/FLAMES campaign are 
displayed as blue filled dots and the position of Lambda Orionis 
(roughly marking the center of the cluster) is highlighted with a 
large blue five-pointed star. 



4.2. Variability connected to activity 

Some hints on the chromospheric activity among the late-type 
population of C69 have already been studied in Paper I through 
the analysis of the alkali lines. Another clear signpost of activity 
in young stellar and substellar objects is variability. This vari- 
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Table 2. Equivalent width (negative for emission, positive for absorption) of the main lines observed in the spectra obtained with 
Magellan/MIKE and Keck/LRIS in its low and high resolution modes (the three different setups are separated by a double line). We 
only show objects where at least one of the forbidden lines, and/or H/3, and/or Hel have been measured. The measurements of Ha 
and the Ca II triplet are shown in Table [6] (with the opposite convention: positive for emission). The symbol "-" indicates that the 
line has not been detected while a white space is indicative of the line falling outside of the wavelength coverage of the spectra. 
Spectral types (SpT) from Paper I (averaged when more than one spectral type was given). 



UDjeCl 


1 




[U1J 


LU1J 


tuij 


TMTTl 


[1M1J 




|M1J 


r^TTi 
LM1J 








SS77 


OJUU 




OJ+O 


OJol 


DO / o 


D / 1 / 


O / J 1 






EW(A) eEW 


EW(A) eEW 


EW(A) eEW 


EW(A) eEW 


EW(A) eEW 


EW(A) eEW 


EW(A) eEW 


EW(A) eEW 


EW(A) eEW 


LOri050 


M4.0 


-9.41±1.64 


-0.53±0.04 


-0.58±0.03 








-0.44±0.07 






LOri055 


M4.0 


-5.11±0.3O 


















LOri056 


M4.5 


-5.16±0.78 


















LOri06() 


M4.0 


-4.13±0.47 




-1.14±0.42 


0.91±0.10 












LOri063 


M4.0 










0.73±0.05 


-2.51±0.39 




-0.93±0.22 


-0.78±0.07 


LOri068 


M4.75 


-5.49±0.62 


















LOri075 


M5.25 




-5.47±1.53 


-1.76±0.50 














LOri098 


M5.0 














-2.24±0.58 






LOril07 


M6.0 






- 


-0.68±0.13 












LOrill4 


M6.0 






-2.39±0.86 






— 








LOrill5 


M5.0 






-0.95±0.03 














LOril39 


M5.75 






-2.76±0.28 








-0.88±0.18 




-1.13±0.30 


LOril40 


M7.0 






-11.82+3.64 








-4.24±0.17 






LOril55 


M8.0 










-4.28±0.82 


-4.42±0.55 


-2.81±0.64 






LOri068 


M4.75 














-1.42±0.24 


-0.76±0.38 


-0.41±0.18 


LOri071 


M5.0 










-0.49±0.11 


-1.13+0.59 






-0.07±0.01 


LOri077 


M5.0 










-0.72±0.25 






-1.14±0.52 




LOri089 


M5.0 










-0.42±0.11 






-0.82±0.38 




LOri091 


M4.75 












-0.11+0.01 


-0.95±0.31 


-0.95±0.69 


-0.34±0.10 


LOri094 


M5.5 














-1.87±0.39 




-1.21±0.82 


LOri099 


M5.5 










-0.58±0.12 










LOril06 


M5.5 










-0.58±0.19 




-1.01 ±0.25 


-0.19±0.02 


-0.26±0.04 


LOril09 


M5.75 














-1.14±0.22 




-0.15±0.04 


LOrilB 


M5.5 












-0.17±0.05 


-1.10±0.14 


-1.73±0.42 




LOrill9 


M5.5 










-4.16+1.03 




-0.58±0.19 


-0.25±0.08 




LOril29 


M6.0 














-0.68±0.15 







ability has been observed not only in the continuum but also in 
a variety of lines, and its dependence with the spectral type has 
also been addressed in the literature (examples can be found in 



Soderblom et al . 1993 Staufferet al. 1997 Barrado y Navascues 
200 1 1 > . In Paper I, we reported a large fraction (>35%) 



et al. 



of the members of C69 showing such variability in alkali lines 
(lithium and sodium). 

Most of the sources showing variability in the alkali lines 
also show variability in the Ha emission, and in Paper I we al- 
ready demonstrated that the differences we see in, for example 
the EW of the lines, are larger than those we would expect just 
because we are comparing observations with different spectral 
resolutions. 

On the other hand, we find a couple of objects exhibiting 
significant variability in Ha and no variability in the alkali lines 
(see Table[4j the top set of sources). The most extreme example 
is LOri073; this object not only shows variations in the mea- 
sured EW of Ha, but also in the profile of the line. As can be 
seen in Fig.[4j while we observed a single peak line profile with 
VLT/FLAMES, a double peak structure was observed several 
years before with Magellan/MIKE. This more complex profile is 



very similar, although narrower, to that presented in Fernandez 
|et al.| ( |2004| l for the weak-line T Tauri star V410 Tau: we can see 
the narrow emission peak slightly red-shifted from the rest ve- 
locity, as well as the shallower component and the blue-shifted 
absorption that could suggest the presence of a wind. 

Finally, another three particularly interesting variable objects 
are discussed in Appendix A. 2.1 They present the peculiarity 
that according to one of the Ha measurements, the emission 



L0ri073 



3 - 




3 - 



Velocity (km/s) 

Fig. 4. Ha profile of LOri073 observe at high resolution in two 
different epochs: with VLT/FLAMES (in red, R-8600) in Jan 
2008 and with Magellan/MIKE (in black, R~ 11250) in Dec 
2002. 
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Table 3. Equivalent width of the main lines observed in the spec- 
tra obtained with VLT/FLAMES. We follow the same conven- 
tions as in Table [2] 



Object 


SpT 


[Nil] 


Hel 


[SII] 


[SII] 






6581 


6678 


6717 


6731 






EW(A) eEW EW(A) eEW EW(A) eEW EW(A) eEW 


C69-IRAC-006 M3.5 


^Q-i-O 0£ 
-u.oy±u.uo 




■If. jVili. 1 j 


-U.J 1 ±U.U J 




M2 5 




-U.4 o±U. 1 J 




n 90-1-0 f\A 
-U.ZU±U.U4 




M3.0 




n a ^ -i- o 1 1 

-U.4j±U. 1 J 






ijWt 1U i J 


M4.0 






-u.zo±u.u / 




T Ori04S 


M3 


n /lO-i-O O/l 








T OriO^O 

ijwl hj JW 


M4.0 


-U.<4O±U.U0 


-U.Oj±U.Z'4 






LOri()53 


M5.0 






9S-1-0 0^ 
-U.ZJ±U.Uj 




LOri059 


M4.5 


-U.0Z±U. 1j 


-KJ.J l±U.Uj 


-U.Oo±U. 1 J 


"3^-1-0 OS 
-u.jO1U.U0 


LOri069 


M5.5 


-0.76±0.13 




-0.46±0.11 


-0.50±0.10 


LOri073 


M5.25 


-0.37±0.09 


-0.22±0.03 






LOri075 


M5.25 


-0.73±0.17 








LOri079 


M6.25 


-1.43±0.07 


-0.40±0.06 




-0.67±0.04 


LOri086 


M5.0 


-1.12±0.13 


-0.34±0.05 






LOri087 


M5.0 


-0.58±0.13 


-0.51±0.03 


-0.90±0.05 


-0.41 ±0.08 


LOri088 


M5.5 


-1.64+0.13 


-0.51±0.09 


-1.59±0.37 


-0.72±0.18 


LOri091 


M4.75 


-0.44±0.06 




-0.81±0.36 




LOri093 


M5.5 


-1.00+0.06 


-0.24±0.03 






LOri094 


M5.5 


-1.71±0.21 


-0.41±0.07 




-0.98±0.19 


LOril05 


M6.0 


-2.68±0.52 






-0.72±0.19 


LOril09 


M5.75 


-2.64±0.44 




-1.42±0.17 


-1.04±0.13 


L0rill2 


M6.5 


-1.76±0.20 




-2.07±0.25 




L0rill5 


M5.0 




-0.83±0.08 


-0.98±0.29 




LOril20 


M5.5 


-3.31+0.17 




-4.73±0.55 


1.34+0.12 



is too large to have a purely chromospheric activity origin (see 
Table|4 second sub-set of objects), but for all the other measure- 
ments this is not the case. 



4.3. Ha emission as a proxy for accretion 

We note once again that the Ha emission can have different ori- 
gins, and although it is commonly used as a proxy for active 
accretion, some considerations have to be taken into account de- 
pending on the nature of the source and its surroundings. 

First of all, and since cool objects are known to be very 
active, the Ha emission can be chromospheric and not related 
to accretion processes. Some limits to this chromospheric con- 
tribution have been suggested in the literature: White & Basri 
42003} proposed EWs of 10A and 20 A for spectral ty pes K3-M2 
and later than M2, respectively. On the other hand, Bar rado y| 
Navascues & Martin ( |2003] > proposed a more spectral depen- 
dent relationship mainly focused on late K, M and L dwarfs 
(and therefore more suitable for our study). This empirical cri- 
terion is based on the saturation limit of chromospheric activity 
([L(Ha)/Lboi] = -3.3) and it is shown in Fig|5]as the boundary 
to discriminate between accreting and non-accreting stars and 
brown dwarfs. 

On the other hand, when trying to extract accurate measure- 
ments from the Ha emission to disentangle between activity and 
accretion, we are also forced to take into account the environ- 
ment where the sources are located. As mentioned before, clus- 
ters like C69, possess a non-negligible nebular component (see 
Pig. [3j, and therefore one has to make sure that the nebular com- 
ponent is subtracted properly from the spectrum of each science 
target. This is not an issue, for instance, when dealing with long- 
slit spectroscopy (like most of our campaigns) but can induce 
some bias when observing with fiber spectrographs if the "sky 
fibers" are not placed carefully. 



In Fig. |3"1 we show the contours corresponding to an Ha im- 
age (from Finkbeiner 2003, with low spatial resolution ~6') of 
the LOSFR, where we have highlighted with crosses the sources 
observed with VLT/FLAMES and with blue filled circles the lo- 
cations of the sky fibers used to correct our measurements. At 
first sight, some structure can be inferred in the nebular emis- 
sion surrounding our science targets. To better characterize this 
effect, we studied the variations of the Ha nebular emission with 
the sky fibers: we measured a mean full width at 10% of ~41 
km/s with a standard deviation of ~3 km/s. We computed a me- 
dian sky fiber and corrected the science spectra with that median; 
therefore the dispersion measured in those fibers translated into 
an added ~7% uncertainty in our measurements. 



The original accretion criterion provided in Barrado y 
Navascues & Martin ( |2003[ ) shows the limiting EW(Ha) as a 
function of the spectral type. Since for most of the sources from 
Dolan & M athieu ( |1999| |2001| l and a fraction of those from 
Sacc o et al.||2008| we do not have spectral type determination, 
we have used the temperature scale derived in Paper I to trans- 
late the original criterion into a T e g vs. EW(Ha) relationship. 

As mentioned in Section |4. 1[ from the total census of 
172 spectroscopically confirmed members of C69, there are 16 
sources for which we do not have a measurement of the equiv- 
alent width of Ha and therefore we cannot apply the criterion. 
Those 16 sources present very dispersed properties such as ef- 
fective temperature, disk presence, etc; and therefore they should 
not produce any bias in the statistics derived for the whole clus- 
ter. 

In Fig [5] we show the accretion criterion applied to the 156 
confirmed members with measurements of the EW(Ha). In or- 
der to be consistent; for objects having a spectral type deter- 
mination, we have also translated them into effective tempera- 
tures using the temperature scale from Paper I (as we did for 
the criterion itself). On the other hand, for objects without spec- 
tral type determination, we have assumed the T e ff derived from 
the SED fit performed with VOSA (|Bayo et al.|2"008] 2012, sub- 



mitted). Applying the saturation criterion, we classified 9%% of 
the members as accretors (14 red dots in the figure; see column 
"Acc" from Table |6]). 

To relate this percentage to the disk presence, in Fig [5] we 
have also included the information regarding the infrared class 
(as a proxy for the presence of disk, from Barrado y Navascues 
et al. 2007 and Morales-Calderon et al. 2012, in prep.) with 



larger circumferences surrounding Class II sources and/or those 
with a MIR slope compatible with thick, thin or transition disks. 
If we only consider the 37 members showing signposts of har- 
boring disks and with available measurements of EW(Ha), we 
estimate that 38^ % show active accretion according to the sat- 
uration criterion. 

To focus now on the substellar population of C69; if we con- 
sider substellar those sources with estimated masses lower than 
0.1M Q to take into account uncertainties in the mass determina- 
tion, we find eight brown dwarfs harboring disks and three of 
them to be accreting according to the criterion. This leaves us 
with a fraction of 37.5+}? % of accretors among the disk-beating 
BDs, very similar to that of the stellar population and compara- 



ble with the substellar one provided in |Scholz et al. ( 2007 1 for 
the similar age cluster Upper Sco (31%, 4 our of 13 objects). 

In Table [6] we present the measured EWs of the Ha line 
of the sources in the sample. For those classified as accretors, 
we have estimated the mass accretion rate using the measured 
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Table 4. Parameters measured and calculated for those sources showing variations in Ha and no significant variability in alkali 
lines; or variations in Ha that causes them to oscillate between accreting and non-accreting (see text for details). Positive equivalent 
widths for Ha denote emission while for the alkali lines it means absorption. We compile also the values from Sacco et al. (2008 ) 
regarding Ha and Li I. 



Name 


EW(Hff) 


EW(Ha) S08 


Ins. Res 


Mass(M ) 


Class 


Disk 


Li I 


Li I S08 


KI 


Nal 


LOri073 


-10.91±1.04 
-4.54 ±0.19 




MIKE R~ 11250 
FLAMES R-8600 


0.34 


III 


Diskless 


0.62±0.03 
0.66±0.05 






2.5±0.02" 


DM048 


-3.36+0.1 1 
-7.72 ±0.19 




CAFOS R-600 
TWIN R-1100 


0.25 


II 


Thick 






0.29±0.07 
0.49±0.21 


2.57±0.43 
2.36±0.18 


LOri068 


-6.82±0.52 
-14.39±0.74 


-7.9±1.26 


MIKE R-I1250 
LRIS R-2650 


0.35 


III 


Diskless 


0.69±0.09 
0.53±0.O5 


0.698±0.025 






LOri075 


-12.81±0.86 
-10.20±0.89 


-10.71±0.98 


MIKE R~ 11250 
LRIS R-950 


0.25 


III 


Diskless 


0.94±0.05 


0.455±0.044 




2.82±0.38 


LOri080 


-21.3±1.00 
-22.82±0.62 
-11.93±1.12 


-14.47±1.25 


BC R-2600 
TWIN R-1100 
LRIS R-2650 


0.15 


III 


Diskless 


0.23±0.07 


0.523±0.020 


1.66±0.89 
1.77±0.33 
2.17±0.46 


2.34±0.46 


LOri091* 


-23.2±12.3 
-13.47+1.11 




TWIN R-1100 
FLAMES R-8600 


0.13 


III 


Diskless 


0.76±0.06 






2.05±0.70 


LOril09 


-19.62±0.52 
-8.70±0.45 




FLAMES R-8600 
LRIS R-2650 


0.16 


III 


Diskless 


0.57±0.11 
0.50±0.15 




2.62±0.57 





"The signal to noise ratio of the TWIN spectrum of LOri091 is very low. This is clear in the very large error on the determination of the EW. Therefore this object is 

not considered variable. 

"Value from |Maxted et al.|(2008) . 



FW, 



et al. 



, %(Ha) and the following relationship derived by |Natta 
d2004f |r 



log(M flcc ) = -12.89(+0.3) + 9.7(+0.7) x 10" 3 FW 10 %(//a) (1) 

The accretion rates calculated in this manner are also listed in 
Table [6] We note that for DM006, classified as accretor, we did 
not have a measurement of the FWk)% and therefore we could 
not estimate the accretion rate. 

In Fig. [6] we show an accretion rate versus mass diagram 
where we compare the values estimated in this work and those 
provided in Sacco et al.| (|2008[ ) also for C69; with those de- 
rived in Rig ia co et al.| ( |201 1[ ) for the slightly younger cluster 
cr Orionis (~3 Myr according to their HR diagram). 

Rigliaco et al. ( |201 1[ ) suggested in their work that there are 
two trends on this diagram with an inflection point at ~O.45M . 
The mass range in C69 for which we have detected active accret- 
ing sources does not allow us to check this feature. Besides, the 
C69 sample is smaller and the dispersion of our measurements 



is much larger than that derived for cr Orionis by Rigliaco et al. 
pOTT) . 

This larger dispersion can arise at a first stage from the differ- 
ent methodology used to estimate the mass accretion rate. While 
Sacco et al. (2008 ) and this work use the Ha emission, "contam- 
inated" by chromospheric activity as already discussed, Rigliaco 
et al. ( 201 1] > use U-band photometry; a "cleaner" methodology. 



On the other hand, C69 is likely older than cr Orionis and the 
accretion disks of the former may be in a different evolutionary 
state than those of the latter. 

Overall, our measurements are consistent with the idea that 
the accretion rate scales with the mass of the central object for 
low-mass stars. But, given the dispersion obtained, this is just 
a very rough trend. Individual sources from Figs. [5] and [6] are 



analyzed in Appendix A. 2. 2 



Finally, to have a better understanding of the relation of the 
Ha emission with the accretion process, we have tried to corre- 



2 The FWio%(Ho-) of the sources are typically above the threshold of 
200 km/s determined by |Natta et al.| (|2004) and both quantities, accre- 
tion rates and FWio%(Lto), are provided in Table|6] 



late that emission with disk properties derived mainly from the 
mid-infrared photometry. 

The theoretical disk models used to interpret the IRAC [3.6]- 
[4.5] vs [5.8]-[8.0] color-color diagram by | Allen et H] p004) 
suggest that the accretion rates increase from the bottom-left to 
the top-right of the Class II region, due to the increase of both 
the disk emission and the wall emission. In Fig. [7] we show the 
mentioned IRAC color-color diagram for the members of C69 
(spectroscopically confirmed members compilation from Paper 
I). We have included information regarding the presence of disks 
(large red circles), the intensity of the emission of Ha (sized blue 
squares) and the classification as accretors (in red). The general 
trend agrees with the disk theory since objects with larger Ha 
equivalent widths (up to accreting sources) have redder colors. 
A similar trend can be observed in the right panel of the same 
figure where we display the mid infrared SED slope as a function 
of effective temperature. Objects with optically thick disks seem 
to exhibit more intense Ha emission. On the other hand, as it was 
already clear in Fig. [5] a large fraction of the sources (~65%) 
harboring disks in C69 do not seem to be accreting from then- 



disks. We will analyze this fact in more detail in section 4.4 



4.4. Disk with low Ha. Binaries clearing the inner disks? 

In the previous section we highlighted the presence of a large 
population of disk-bearing sources that do not show any sign- 
post of active accretion. In Fig. [8] we show examples of SEDs 
corresponding to this class of sources. We looked for some char- 
acteristic that would differentiate these objects from the others in 
our sample (other than the measured Ha equivalent width). Their 
effective temperatures are mainly colder than ~3750K, but, as 
explained in the next section, that is characteristic of all the disk 
population of C69. 

One possibility for these "quiet" disks would be that they 
are dissipating their inner disks, but we see no difference in the 
near-infrared colors with respect to the actively accreting disks. 

Besides, adopting the characterization of the a parameter 
from Lad a et aX] (|2006 ), these objects harbor mainly optically 
thick disks (60% of them), but there are also sources with op- 
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Fig. 7. Confirmed members with Ha in emission compiled in Paper I. Left: IRAC color-color diagram including information about 
the nature of the disk and the Ha emission of the sources. Sources showing infrared excess have been surrounded by a black 
circumference. A blue square with a size proportional to the intensity of the Ha emission of the sources is also show n in the figure. 
Whenever no square is present, the EW of the line is lower than 5 A. The Class II area after Allen et al. (2004i is highlighted 



with dashed lines. Right: Effective temperature vs. IRAC SED slope of the confirmed and candidate members of C69 including 
information about the nature of the disk and the Ha emission of the sources as in the left panel. 



tically thin disks and the so called transition/cold disks (Merin 
|et al.|2010| ). There are several mechanisms to explain the evolu- 
tionary status of these transition disks; one of the most attractive, 
in the context of planet formation, is the clearing of the inner disk 
by a giant planet in its earliest stages of formation (see the first 
observational candidate for this scenario in|Huelamo et al.| 201 



reland 



Another possibility is tidal truncation in close binaries 
& Kraus 2008). This seems to be the mechanism at work in at 
least one of the cases of the low Ha transition disks of C69; 



LOri043 that was classified as SB2 by |Maxted etaT1p008| ). On 
the other hand, LOri043, is the only documented spectroscopic 
binary in this sample of quiet disks. 

Finally, only two of the sources from this set have been de- 
tected in X-rays, suggesting that these are not particularly active 
objects either. In conclusion, we could not find any parameter 
(other than the Ha emission) that unites these objects or differ- 
entiates them from the actively accreting population. 



5. Disk vs diskless populations 

By the end of the previous section, we showed that we cannot 
trace differences between the actively accreting - quiet disk pop- 
ulation of C69. In this section we will analyze the clear distinc- 
tion between disk and diskless sources, in terms of their X-ray 
emission and their mass functions. 



5.1. Ha and X-ray 



As stated in Section 2; Barrado et al. (201 1 1 presented the anal 



ysis of the XMM-Newton observations of two fields in C69. 
Several months later, Franciosini & Sacco ( 201 1[ ) complemented 
the study by adding an extra field that covers the vicinity of the 
massive star A Ori, roughly at the center of the cluster. 



These X-ray observations should trace well the weak-line T 
Tauri (and substellar analogs) population of C69; therefore we 
have combined the data from the two studies and correlate it 
with our census of spectroscopically confirmed members. 

In Fig. [9] we illustrate this advantage of the X-ray observa- 
tions to unveil the weak-line T Tauri population. We show every 
member of C69 lying in the field covered by XMM-Newton ob- 
servations. We have highlighted in blue objects above the com- 
pleteness limit of 0.3M Q . We see how most of the sources clas- 
sified as weak-line T Tauri, according to the saturation criterion, 
are detected in X-rays. We also show that the objects that are not 
detected in X-rays are preferentially those harboring optically 
thick disks (both, active accretors and non-accreting sources). A 
total of 9 sources out of 12 with masses above O.3M0, harbor- 
ing disks and within the XMM-Newton fields of view are not 
detected in X-rays and are labeled in Fig. [9] 



This dichotomy is shown even clearer in Fig. 10 where we 
show the X-ray luminosity vs bolometric luminosity ratio as a 
function of Ha equivalent width. Here, we see how most ob- 
jects with EW(Ha) between ~5-20 A and not harboring disks 
are detected, while those showing infrared excess are not. We 
have highlighted five of the non-detected sources discussed in 
the previous paragraph for which Barrado et al. (201 1 ) provide 
upper-limits of the X-ray luminosity. The remaining four non- 
detections are located towards the center of the cluster, within the 
field of view of the observations by Franciosini & Sacco ( 201 1) 
and no upper-limit for the X-ray luminosity of the sources is 
provided in that work. 

5.2. Disk fractions 

Regarding the formation mechanism of brown dwarfs; if they are 
just a scaled down version of stars, one would think that the disk 
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Fig. 9. Spectroscopically confirmed members from Paper I that lie in one of the two XMM-Newton pointings from |Barrado et al. 
( |201 l| l or in the central one from [ Franciosini & Sacco|(|201 1[). We show in different colors objects with mass below (in black) and 
above (in blue) 0.3MO (the completeness limit for Barr ado et al.||20lT| ). For the EW(Ha) measurements we have gathered data 
from Sacco et al. ( 2008|l; Dolan & Mathieu ( 1999]> and this work. We have highlighted with large stars those sources detected in 
X-rays according to Barrado et al.| ( |2011| l and/or |Franciosini & Sacco| ( |201 1 1, and we have joined with a grey solid line different 
measurements of EW(Ha) corresponding to the same source. Left: Effective temperature vs Ha equivalent width (positive for 
emission) along with the saturation criterion. Right: Mid-infrared slope vs Ha equivalent width (same convention as in the previous 
panel) 



life-times above and below the hydrogen burning limit should be 
the same. On the other hand, there is some evidence for longer 
life-times for infrared excesses in very low-mass stars and brown 
dwarfs than in higher-mass stars (Lad a et al.|2006| Allers et al. 
|2007) . 

To test these possible differences, instead of estimating a typ- 
ical star and brown dwarf disk fraction, we have divided our cen- 
sus of confirmed members into two sets: sources showing some 
infrared excess, and objects detected also in the mid-infrared, 
but exhibiting purely photospheric colors. The first set contains 
43 sources and the second one 115. We propose this approach to 
avoid choosing a fixed frontier in mass between stars and brown 
dwarfs so that our findings are easier to compare with other stud- 
ies. Regarding the sources with excess, we have considered, in 
the same set, all kinds of MIR slopes and shapes: optically thick, 
optically thin and transition disks. 

We have computed the cumulative mass functions for both 
samples and the result is shown in Fig.[TT] there is a clear differ- 
ence at ~0.6M Q (~M2 spectral type). While the mass function of 
the diskless population rises up to ~1.1M , the disk population 
is almost completely composed of sources with masses lower or 
equal than O.6M . In other words: while we find diskless objects 
for every bin in mass, sources more massive than ~O.6M seem 
to have lost their disks already. We tested the dependence of this 
change of behavior with the 5 Myr age assumed for C69 to esti- 
mate the masses. While an older age will significantly affect the 
mass determination of the lowest mass members of C69; that is 
not the case for sources with masses above 0.3M Q according to 
the isochrones by |Baraffe etal] ( [T998) ; |Chabrier et IT] ([2000 ); 



Allard et al. ( 2003) , and therefore our result is robust against 
changes in the age determination of C69. 

On the other hand, if we use these two cumulative fractions 
but without normalization to estimate the disk fraction as a func- 
tion of the mass, we can see that the situation is more complex. 
In Fig. [12] and Table [3] we show the ratio of the two cumulative 
functions; that is, for a given mass M,, we provide n,-/iV,- where 
n, is the total number of sources with mass < M,- and infrared 
excess, and Nj is the total number of sources with mass < M, 
(rij plus the sources with purely photospheric infrared slopes and 
masses < Mi). For the error treatment, since our sample is very 
large, for masses larger or equal than 0.3M Q we can derive stan- 
dard Poisson uncertainty limits and for the lower masses we have 
used the approach described in Burgasser et al. ( 2003) . 

Regarding the completeness of the diskless population, even 
though the X-ray observations are only complete down to 0.3M G 
we are confident that this does not induce any bias in our analy- 
sis. The reason is that every source that in the spectroscopic con- 
firmation turned out to be later than M3.5 and was detected in 
X-rays, had previously already been selected as an optical pho- 
tometric candidate. Therefore the optical photometric selection 
was as good as the X-rays one picking up the low to very low- 
mass members of C69. 

The disk fraction function seems to peak (77%) at the brown 
dwarf boundary, dropping abruptly with increasing mass up to 
~0.3M o . It then stabilizes at ~33% before falling again for 
masses higher than ~0.6M G . The total stellar disk fraction is 

26 + _\ % 
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Fig. 5. Ha equivalent width versus the effective temperature 
for confirmed members of C69. We display data from 
& Mathieu] (|1999[ [200T) with filled squares, and from 



Dolan 



Sacco 



et al.| ( |2008| ) with filled diamonds. Our data are displayed with 
solid circles. In every sample, red symbols are used for those 
sources classified as accretors. Overlapping large circles high- 
light sources exhibiting excess in the Spitzer/IRAC photometry. 
For some sources we had more than one epoch of data (either 
ours or from |Dolan & Mathieu|1999l|200T] |Sacco et al.|2 008 >. 
These different measurements for individual objects appear joint 
with a solid line. The short-dashed line corresponds to the satura- 
tion criterion defined by Barrado y Navascues & Martin (2 003) >. 
A vertical dotted line highlights the substellar frontier for an es- 
timated age of 5 Myr according to the isochrones from the Lyon 
group. Particular sources discussed in Appendix A.2.2 are high- 
lighted with grey labels. 



We must note that the extremely high disk fraction for sub- 
stellar objects should be taken as an upper-limit. In this low mass 
regime, we have some sources that have not been detected in the 
two reddest channels of Spitzer/IRAC (5.8 and 8.0 fim). We did 
not consider those sources since we only used sources with pho- 
tospheric MIR SED to estimate the disk fractions. If we were 
certain that those sources do not show excess (and the fact that 
they are not detected at those red wavelengths is a good indicator 
of that), the percentage would decrease down to ~58% 



This fraction is still larger than that derived by Barrado y 
|Navascues et a l. (2007 ) for the same cluster (~40%). This is not 
too surprising since in this study we consider the fraction with 
any kind of disk and not only Class II sources as in Barrado y 
|Navascues et al.| ( |2007|l. On the o ther hand, the difference with 
the value derived by jScholz et al.| ( [2007] i for Upper Sco (-37%), 
a similar age cluster; although significant, could just be caused 
by small number statistics or again, use of different criterion to 
infer whether a source is harboring a disk or not. 
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Fig. 6. Mass accretion rates versus mass of the central object for 
accreting objects. In red and blue we show members of C69 (red 
dots for measurements from this work and blue squares for those 
from |Sacco et al.|2008| the masses of LOri061 and LOri063 have 
been slightly shifted so that in the figure the comparison between 
the measurements of two studies is clearer) and in black mem- 
bers of the slightly younger cluster cr Ori (measurements from 
Rigliaco etal.|20lT] ). We have included blue and red vertical bars 
to highlight special cases where more than one measurement is 
available for a given source (see Appendix |A.2.2| for details). 
In all cases, the masses shown in this figure are the average be- 
tween the one derived from the Lboi and the one derived from the 
T e ff. The error bars display the differences among these determi- 
nations, see Paper I for details. Particular sources discussed in 
Appendix A.2.2 are highlighted with grey labels. 



Table 5. Disk fraction as a function of mass. 
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5.3. Spatial distribution 

In Fig.[T3]we show the spatial distribution of C69 spectroscop- 
ically confirmed members (by |Dolan & Mathieu| 1999| [2001 ; 
Sacco et al.|2008|[Maxted et al.|2008|or this work) including in- 



formation about X-ray emission (from Barrado et al. 201 1 and/or 
|Franciosini & Sacco 2011), the presence of a disk ( |Barrado y| 



[Navascues et al.|2007 1 and whether the disk is accreting or not. 
As can be seen in the figure, the sources with a disk show a 
higher concentration towards the center of the cluster with re- 
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Fig. 8. SEDs of the sources classified as Class II based on their 
IRAC colors and showing EW(Ho-) low enough not to be clas- 



sified as accretors according to Barrado y Navascues & Martin 
< |2Q03> . 
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Fig. 10. EW(Ha) vs LxA-boi diagram for the members from 



Paper I in the area covered in X-rays by either |Barrado et al. 
(201 1 1 or Franciosini & Sacco (2011 1. As in previous figures we 



use large open circles to highlight objects showing mid-infrared 
excess and red dots for those classified as accretors. Five pointed 
stars indicate that the estimated mass of the object is larger than 
the critical mass 0.6M G explained in section[5^2| 



Fig. 11. Cumulative (normalized) mass functions for sources 
showing signs of harboring disks (infrared excess), plotted with 
large circles surrounding the points of the histogram, together 
with the cumulative mass function of sources detected by IRAC 
showing photospheric MIR slope. The vertical red dashed lines 
highlights the mass (~O.6M ) at which the disk population is 
almost complete. 
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Fig. 12. Disk fraction as a function of the mass of the members 
of C69. For the ratio we use the same sets of objects as in the 
previous figure. 
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spect to the diskless population (contrary to what one would 
expect according to the supernovae scenario and already sug- 
gested in Barrado y Navascues et al. 2007| ). In fact, if we assume 
that actively accreting systems are younger than those which do 
not show any sign of accretion, the youngest population of C69 
seems to be clustered either around the central star A Ori or to 
the Soufh-West. 

We have computed the two-sided Kuiper statistic (invariant 
Kolmogorov-Smirnov), and the associated probability that any 
of the previously mentioned populations were drawn from the 
same distribution. The tests reveal that the cumulative distribu- 
tion function of Class II candidates is very different from that of 
Class III, with a 99.9% probability that both populations have 
been drawn from different distributions. 

The conclusion from this test is that objects with masses 
lower than O.6M have been less efficient, in the life time of 
C69, in loosing their circumstellar material. 
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Fig. 13. Spatial distribution of the spectroscopically confirmed 
members of C69. Active accretors according to the saturation 
criterion are highlighted in red. Information regarding infrared 
excess (as a proxy for the presence of disk, from Barrado y 
Navascues et al.|2007 1 is provided by surrounding those sources 
with larger circumferences. The more massive population of C69 
(A Ori itself and the B stars) is shown with grey four-pointed 
stars. Finally, the members with X-ray detections are shown with 
five-pointed blue stars. 



6. Summary and conclusions 

We have analyzed the spectroscopic properties of the very com- 
plete sample of members of C69 compiled in Bayo et al. ( 201 l| l. 
Using different spectroscopic features we have tried to better un- 
derstand the similarities and differences between the disk and 
diskless populations and the stellar and substellar ones. Our main 
results can be summarized as follows: 



We have estimated the rotational velocities for eight mem- 
bers. We find a high dispersion in the vsin(/) values, being 
larger among the diskless population. We interpret this as a 
result of disk locking in some of the C69 members. 
We have studied the variability of the Ha emission line 
in 142 members classified as magnetically active (non- 
accreting sources). We have also identified candidates that 
might have experienced flares during the epoch of observa- 
tion. 

We have tried to disentangle the activity and accretion contri- 
butions to the Ho- emission. In this context, for those sources 
showing large Ha equivalent widths (larger than the one ex- 
pected to arise from chromospheric activity) we have derived 
accretion rates using mainly the full width at 10% of the flux 
of the Ha emission line (but also the Ca II IRT) obtaining a 
very large spread of values. Once individual cases are ana- 
lyzed, these spread values are still compatible with those pre- 
viously reported in the literature for objects of similar mass 
and age (Riglia co et al.|20lT l. When considering all the con- 



firmed members with Ha measurements, we estimate a 9*1 
% fraction of accretors in C69. If we only consider objects 
with disks, 38^ % show active accretion. 
We have studied the relation Ha - disk properties. While the 
general trend expected from disk models applies, we also 
identify a pretty large population of "quiet" disks. Objects 
showing clear mid-infrared excess (with a variety of SED 
shapes) but Ha levels compatible with arising from pure 
chromospheric activity and not from the interaction with the 
disk via accretion. 

We find a stellar disk fraction of 26^ % and we can put some 
limits to the substellar disk fraction (for the faintest mem- 
bers, due to sensitivity limitations, the lack of detection at 5.8 
and/or 8 micron does not imply the presence or absence of a 
substellar disk) of 58%. These fractions do not compare too 
badly (taking into account the difference in procedure fol- 
lowed to estimate the fraction) with that previously provided 



for C69 itself by Barrado y Navascues et al. 



and is significantly higher than that derived by Scholz et al. 



(2004b) (-40%) 



(2007 ) for the cluster of similar age Upper Sco (~37%). 
Regarding the accretion fraction in the substellar domain, we 
do not see dramatic changes from the global fraction of C69. 
This accretion fraction for brown dwarfs varies from 30% 
to ~43% according to our uncertainties in mass determina- 
tion, This range is also compatible with the one provided by 
|Scholz et al.| ([2007 ) for the substellar population of Upper 
Sco. 



We have confirmed that X-ray observations (Barrado et al. 
|201 1| |Franciosini & Sacco|2011[ ) are extremely efficient re- 
covering the Class III population of C69 in the intermediate 
to low-mass range. 

We have compared the mass function of the disk and diskless 
populations of C69 finding that O.6M seems to be the criti- 
cal mass below which a significant fraction of the members 
preserve their disks. This result implies different disk life- 
times for different stellar masses. In particular, for masses 
lower than 0.6M Q we have shown that the disk fraction rises 
very steeply with the caveat that in the brown dwarf domain 
the fraction provided should be taken as an upper-limit. 
We have studied the spatial distribution of the disk-harboring 
population of C69 and we find that, opposite to what Dolan 
|& Ma thieu (2001) derived from their more massive mem- 
bers sample; the density of disk-sources is larger closer to 
the center of the cluster, which is inconsistent with the SN 
scenario invoked to explain the origin of the Lambda Orionis 
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star forming region. In addition, the winds (current or in the 
recent past) from the massive star A Ori, seem not to have af- 
fected the distribution of disk and diskless cluster members. 
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Appendix A: Particular sources 

In this appendix we provide further analysis for sources show- 
ing peculiarities in the properties studied in Sections [3]and[4] To 
keep consistency throughout the paper we have grouped the in- 
teresting sources from Section [4] following the same subsection 
scheme. 



A.1. Rotational velocities 

LOri075: This source has been classified as single-line spec- 
troscopic binary (SB1) by |Maxted et al. (2008) (but no binarity 
signpost has been reported in |Sacco et al.| |2008). According to 
|Maxted et al.| ( |2008| l, the spectral lines for this star show rota- 
tional broadening; they compared them to those of a narrow- 
lined star of similar spectral type and estimated a projected rota- 
tional velocity of vsin(/) ~ 65 km/s. They classified the source as 
SB1, but they also noted that there is an asymmetry in the cross- 
correlation function (CCF) in the form of a blue-wing, partic- 
ularly when the measured radial velocity corresponds to a red- 
shift. Therefore they suggested that the fainter component in this 
binary was detected but unresolved in their spectra. 

We have detected a double peaked structure in Ha and Li I in 
our Magellan/MIKE spectra which made us believe that we had 
spectroscopically resolved the source. While the origin of double 
peak in Li I should be related to binarity, the Ha one could arise 
from an accreting companion for example. Further research on 
the structure of some photospheric lines marginally confirms this 
idea. In Fig A. 1 we show the double peaked structure found in 
some of the most prominent photospheric lines for this object 
(given the low temperature of the source, these "most prominent 
lines" are still very weak). We have measured a mean relative 
velocity of ~45 km/s (<x ~ 9 km/s). We have synthesized a 3500 
K (log(g) = 4.0) Kurucz spectra (an effective temperature ~100 
K hotter than the one estimated for our source, but the coolest 
temperature for the Kurucz collection) in the region of the Ba 
/15535A line with the same resolution and a rotational velocity 
close to the one derived by |Maxted et al. (2008 1 (~50 km/s). We 
have checked that the closest line in the synthetic spectra has a 
relative velocity of ~1 10 km/s, much higher than those measured 
by us. We show on the right-hand side panel that the relative 
velocity derived for the photospheric lines does not agree with 
the one that would be measured from the Ha profile. This fact 
and the weakness of the lines measured force us to consider the 
resolution of the binary as tentative. We must note anyway that 
the environmental Ha component (see Fig [3]) of the region or a 
possible accreting companion could change the relative velocity 
of the peaks of this emission line. 



A.2. Activity and accretion 

A.2.1. Variability connected to activity 

In our study of the Ha variable sources, we find a sub-set of 
five objects for which the criterion from Ba rrado y Navascues &| 
Martin (2003) applied to spectra taken at different epochs pro- 
vides contradictory results. While for some measurement of the 
object taken in one epoch (EW(//a)i), the source would be clas- 
sified as accretor; for a different epoch measurement of the same 
object (EW(Hff)2), the Ha emission could be explained purely 
in terms of activity. 

- LOri068 and LOril09 were observed twice during our 
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Fig.A.l. Left: Double peaked structure found in photospheric 
lines for LOri075. Rest frame velocity and mean relative ve- 
locity of the second peak are indicated with dashed lines. The 
shaded (blue) rectangle shows the ±3cr area of the second peak 
location. Right: Ha profile of the same source. Note how the 
secondary peak dashed line location (calculated from the photo- 
spheric lines) does not agree with the position of the peak (see 
text for details). 



campaigns and LOri068 was also observed by Sacco et al. 
(2008). Both objects are classified as diskless sources based 



on their IRAC slopes, and while our Li I measurements agree 
among themselves (and for the case of LOri068, with the 
one provided by Sacco et al.|2008| l, in one spectrum for each 
source, the Ha emission is much more intense than in the 
others. We believe those spectra were taken while the ob- 
jects were experiencing a flare. In Fig A. 2 we compare the 
line profiles of those objects in "steady" and "flared" states; 
and we show that, while in the case of LOri068, we can see 
some asymmetry in the line profile for the intense emission 
(which would indicate mass motion), that is not the case for 
LOril09 (even though for the latter the change in EW is 
much stronger). 

LOri091 cannot be classified with certainty as a variable 
source. There is one measurement of Ha clearly off from 
the other two available, but that measurement corresponds 
to a TWIN spectrum of very poor S/N which translates into 
a large uncertainty on the continuum, and therefore a very 
large error-bar in the measurement. 

LOri075 is an un-resolved (or marginally resolved, see sub- 
section |AT| double system, and therefore variability in the 
measured Ha is not surprising. 

LOri080 is a puzzling case. We observed the object twice; in 
2003 (at Las Campanas) and 2005 (at Calar Alto; see Paper 
I for a description of the instrumentation used in each case), 
and both measurements agree within the errors. These mea- 
surements place LOri080 at the border of being classified 
as accreting according to the saturation criterion (see next 
subsection), although the object shows no infrared excess 
in the IRAC data. No peculiarity has otherwise been found 
regarding the profiles of the lines in either spectra. On the 
other hand, Sacco et al. (2008 ), also observed LOri080 with 
FLAMES and found an EW for Ha significantly lower. 
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Fig.A.2. Comparisons of the Ha line profiles for different ob- 
servations of the same two objects: Upper panel: LOri068; 
one of the objects we suspect experimented a flare during the 
Keck/LRIS observations (R~2700). Some asymmetry can be 
seen in the line even though the resolution of the spectrum is 
moderate. Lower panel: LOril09; another object suspected to 
have a flare for which no asymmetry in the line has been found. 



A.2.2. Ha emission as a proxy for accretion 

In the process of determining accretion fractions, ratios and their 
relation with the disk properties, we encountered several partic- 
ular cases that we describe below. 

LOril61: Is the brown dwarf from Fig[5]that, even though its 
Ha emission places it well above the saturation criterion, has not 
been classified as harboring a disk according to its SED. The is- 
sue with this very faint source is that it was not detected in IRAC 
channels three and four (5.8 and 8.0 micron, respectively). Since 
the sensitivity of these channels is lower than that of one and two 
(3.6 and 4.5 micron), it could be the case that this object indeed 
has a disk that we are not sensitive to and that is undergoing ac- 
cretion. In that scenario, the estimated accretion rate according 
to the FWio%(fla) would be ~ 1 x 1(T 10 M e /yr, which is much 
lower than the accretion rate derived for LOril56 (~ 9.5 x 10~ 9 
M Q /yr), also a brown dwarf with the same spectral type and dis- 
cussed later on in this section, but harboring an optically thick 
disk. 

An example of such a disk would be a transitional disk, 
where the excess would be only detectable at larger wavelengths. 
We checked the new release of the WISE catalog (in the prelim- 
inary version the source is not detected) and we found a coun- 
terpart within 1". Unfortunately, although the photometry at the 
largest wavelengths (~11 and 20 micron) shows a clear excess, 
these measurements have been classified as "U" (upper limit), 
and therefore we cannot confirm that this source does harbor a 
disk. 

C69XE-009: An X-ray candidate from |Barrado"^taLl ( pOTT] l 
confirmed spectroscopically as C69 member in Paper I. This ob- 



ject is right at the limit of the saturation criterion; based on its 
SED it was classified as a candidate transition disk, but the lin- 
ear fit to the mid-infrared slope is photospheric. Given that it is 
clearly an active object (detected in X-rays), and that the disk 
possibility is based on a very slight excess detected only in one 
infrared band, we assume that the Ha emission has its origin in 
chromospheric activity and not accretion. 

C69-IRAC-005: Is the source from Fig. [6] exhibiting the 
largest accretion rate based on the FW[o%(Ha) (log(M) = 
-5.56 + 0.25) and the Spitzer/IRAC photometry suggests that 
it harbors an optically thick disk. This particular source was ob- 
served with CAFOS in low resolution mode, with a wider wave- 
length coverage than the other instruments used (see Paper I). 
Thus, we have been able to obtain a different estimation for the 
accretion rate based on the equivalent width measurement of the 
components of the Call triplet (at 8498 A, 8542 A and 8662 A). 
This emission could be a sign of chromospheric activity too, as 
in the case of Ha; but the obtained equivalent widths for the 
triplet are too large; placing our measurements in the broad-line 
component of the unresolved line structures that is generally 
related to accretion (see Comeron et al.|[2003| llvlohanty et al. 



2005 ). Furthermore, as in Comeron et al. (2003 1, the Call triplet 



line ratios are very close to 1:1:1 (quite different from the 1:9:5 
expected ratio for optically thin emission). 

We used the following equations to estimate the accretion 
rate from the Call triplet (these equations were derived from the 
accretion line pr ofile study by|Muzerolle et al.|19 98 and are fur- 
ther discussed in Comeron et al. 2003 ): 

log(M flce ) = -34. 15 + 0.89 log(Fcaii(i8542)) (A. 1) 

Fcanu8542) = 4.72xl0 33 £W(Can(^8542))xl0" a4(m ^ a54Al ' ) (A.2) 

where i*caii(,t8542) is the flux in the line, m A is the magnitude of 
the star at /18542, and Ay is the visual extinction translated to 
the wavelength of the line of study using Fitzpatrick ( 1999) rela- 
tions. Since the bluest photometric point that we had for this ob- 
ject is the 2MASS J magnitude, we used the best fitting model to 
the SED of the source as a scaling factor to estimate nix- On the 
other hand, according to the intrinsic colors by Leggett (1992) 
and our determination of the spectral type (M3), we find a very 
low Ay value of 0.03 mag (quite lower than the average value of 
0.36 mag derived for the cluster by |Duerr et al. 1982 but neither 
of them would significantly affect this estimation). 

We obtained an accretion rate value of ~3xl0~ 7 Mo/yr; al- 
most an order of magnitude lower than the one obtained based 
on Ha, which gives us an idea of the caveats of estimating ac- 
cretion rates from measurements that can be well contaminated 
by activity or even by wind contributions. 

On the other hand, even with the two estimations differing 
by such a large factor, this object still seems to be experiencing 
heavy accretion. We have compared its spectrum with that of 
C69-IRAC-002 (another M3 star, observed with the same setup, 
harboring a disk but with Ha compatible with pure chromo- 
spheric activity and no Call emission) looking for veiling emis- 
sion, and no blue excess has been found in the source (further 
than a marginal excess right in the blue edge of the spectra that 
we think corresponds to an instrumental signature rather than a 
real excess). This result is not surprising since the wavelength 
coverage starts at 6200 A, and veiling in young stars is normally 
detected at bluer wavelengths. Therefore, we would need further 
spectroscopic measurements to confirm the presence of veiling 
in this source. 

LOri050, LOri061 and LOri063 are the other sources with 
more than one estimation of M,,,.r. 



A. Bayo et al.: Spectroscopy of Very Low Mass Stars and Brown Dwarfs in the Lambda Orionis Star Forming Region., Online Material p 3 



LOri061 Ha 



C69-IRAC-007 



C69-IRAC-006 



I 0.8 




6520 



6540 



6560 6580 
Lambda (A) 



6600 



Fig.A.3. Detail of the Ho- emission of LOri061 in the 
CAHA/TWIN spectrum. Note the dependence on the measure- 
ment of the full width at 10% of the flux with the pseudo- 
continuum choice, in particular for the light green and the teal 
cases (for a complete description on the process to determine 
the different continuums see Appendix A of Paper I). 



LOri050 is a spectroscopic binary according to Sacco et al. 
(12008b andlMaxted et al.l (I2008JI; it has been classified as Class 



II according to its mid-infrared photometry and we obtained 
two spectra with different instrumentation (see Paper I for de- 
tails). According to the Ha emission, in both cases, the object 
is above the saturation criterion. The estimated accretion rates 
for both measurements agree well within the errors (log(M) = 
-11.09±0.05, -10.91+0.05). Therefore, we are observing a very 
interesting system with a total stellar mass of ~0.3 M Q and a cir- 
cumbinary disk actively accreting. 

For LOri063, on the other hand, the two available accretion rate 



estimations (from |Sacco et al. 2008 and this work) differ by more 
than an order of magnitude (log(M) = -10.7+0.3, -11.87+0.07). 
LOri063 harbors an optically thick disk according to its IRAC 
photometry, and the change in the full width at 10% of the flux 
in Ha is also reflected in the change in EW of the line (>9A). 
Finally, our measurement of LOri061 does not agree at all 
with that from Sacco et al. ( 2Q08) > (two orders of magnitude 
difference, log(M) = -10.2+0.3, -7.65+0.05). We believe this 
difference arises from how sensitive the measurement of the 
FWio%(ffa) is to the local continuum determination. In Fig. A. 3 



we illustrate the case graphically. While our automatic procedure 
(see Appendix A of Paper I for details) identifies a local contin- 
uum, the thick light green line, other routines fitting global con- 
tinuum could base their measurements on the teal line. This dif- 
ference in the determination of the "real base" of the line yields 
the large discrepancy in the estimated accretion rate. We must 
note in any case that among our data on accretors, we do not 
have other sources where the Ha profile can provoke this confu- 
sion in the continuum determination. 

LOril26, LOril40 and LOril56 are the three brown dwarfs 
(LOril26 is right at the limit between BD and very low mass star 
depending on the method used to estimate its mass, see Paper 
I) from Fig. [6] exhibiting very large accretion rates (log(M) = 
-8.78+0.10, -8.88+0.14, -8.02+0.10, respectively). According 
to their mid-infrared slope, the three targets harbor optically 
thick disks. And according to their very large Ha equivalent 
widths, they are well above the saturation criterion. 

With such high accretion rates some veiling (due to excess 
emission from the accretion shock) could be expected in these 
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Fig.A.5. Double peaked structure in the Ha emission of C69- 
IRAC-007 and C69-IRAC-006. 



sources (as it is the case for LS-RCrA 1, Barrado y Navascues 



et al. 2004a I. In order to study this possibility, for each brown 
dwarf, we selected a non-accreting Class III source with very 
similar spectral type (one half subclass) and we compared the 
strength of several TiO molecular bands in both spectra. As can 
be seen in Fig |A.4| no significant differences are found in the con- 
tinuum level of any pair of sources. In fact, in the three cases, r A , 
defined as FWexcess/fWphotosphere is neglig ible. Whilst for LS- 
RCrA 1, Barrado y Navascues et al. (2004a) found that r,\ varies 
from ~1 to ~0.25 for the wavelength range 6200-6750 A we 
find a horizontal slope in this interval. The only cases where a 
linear horizontal r A does not work are located on the very edges 
of the detector, and therefore we can conclude that no veiling is 
detected in any of the spectra in the studied wavelength range 
(this does not imply that some veiling cannot be present at bluer 
wavelengths). 

C69-IRAC-006 and C69-IRAC-007: Both sources are clas- 
sified as Class II based on their IRAC photometry, and they show 
double-peaked structure of the Ha emission as can be seen in 
Fig [A3] Whilst the sky subtraction for C69-IRAC-007 worked 
very well, some residual could remain in the case of C69-IRAC- 
006 (although we do not see any structure on other, very narrow, 
"sky lines"). 

This double peak is not present either in the other emission 
line detected in both spectra (He I) or in the absorption lines, 
discarding to a certain extent the possibility of these sources be- 
ing spectroscopic binaries (SB2). We do not have an estimation 
of the rotational velocities of either given the resolution of the 
VLT/FLAMES observations (R-8000), but these almost sym- 
metrical double peak structures in Ha have been reproduced 
with models for higher mass stars with rapid rotation seen al- 
most pole on (see Muz erolle et al.|20"03 and references therein). 

The accretion rate calculated for C69-IRAC-007 is shown 
in Table [6 ] since this source fulfills the Barrado y Navascues 



& Martmp|2003l l criterion (log(M) = -7.05 + 0.02). On the 
other hand, although the measured Ha equivalent width of C69- 
IRAC-006 lies well below the saturation criterion of Barrado y 
Navascues & Martm| ( |2003| l, the wide FWio% measured (~190 
km/s) places this object very close to the limit of accretors ac- 
cording to Natta et al.|2004 ). In addition, note the resemblance 
of the Ha profile of C69-IRAC-006 with that of Cha Ha2, an 
accreting brown dwarf, modeled in detail (and showing pecu- 
liarities attributed to the presence of an outflow) in Natta et al. 
d2004] >. 
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Fig. A.4. Upper panel: Comparison of accreting (red, very low mass / BD LOril26 and brown dwarfs LOril40 and LOril56) and 
non-accreting (black, with the same spectral type) members with low resolution spectra. Note the absence of veiling. Lower panel: 
Detail of the Ha emission o the same comparison. 
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